The Adiabatic Theorem for Quantum Systems with Spectral Degeneracy 
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By stating the adiabatic theorem of quantum mechanics in a clear and rigorous way, we establish a 
necessary condition and a sufficient condition for its validity, where the latter is obtained employing 
our recently developed adiabatic perturbation theory. Also, we simplify further the sufficient con- 
dition into a useful and simple practical test at the expenses of its mathematical rigor. We present 
results for the most general case of quantum systems, i.e., those with degenerate energy spectra. 
These conditions are of upmost importance to assess the validity of practical implementations of 
non-Abelian braiding and adiabatic quantum computation. To illustrate the degenerate adiabatic 
approximation, and the necessary and sufficient conditions for its validity, we analyze in depth an 
exactly solvable time-dependent degenerate problem. 
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The adiabatic theorem [l[ has played, and still plays, 
a fundamental role in practical quantum physics applica- 
tions. Indeed, the ability to determine how the slow dy- 
namics of external probes coupled to a system affect its 
time evolution has applications ranging from the notion 
of thermal equilibrium and non-equilibrium phenomena 
Q to the conditions under which an adiabatic quantum 
computer can reliably operate Q. Useful and practical 
quantitative conditions for the validity of the adiabatic 
theorem are also relevant to the important current prob- 
lem of assessing the feasibility of any information process- 
ing scheme that uses the concept of fractional exchange 
statistics and non-Abelian braiding 

General physical principles dictate that, in three space 
dimensions, elementary particles can only obey fermionic 
or bosonic statistics. Kinematic constraints do not allow 
for fractional exchange statistics: electrons are spin-1/2 
fcrmions and photons are spin-1 bosons. Nonetheless, 
fractional statistics particles or modes may emerge from 
the collective behavior of elementary particles, i.e., col- 
lective excitations of a quantum field, as a result of a 
dynamical process. The latter requires special circum- 
stances and constraints that should be analyzed on a case 
by case basis. For instance, for two localized degenerate 
Majorana modes to realize a non-Abelian braiding pro- 
cess we need to design the physical Hamiltonians realizing 
the braiding that do not lift the degeneracy and can be 
implemented adiabatically. If those constraints are not 
met experimentally then the braiding operation is faulty. 
Physical systems where such fractional statistics emerges 
have highly degenerate energy spectrum, thus justifying 
a careful statement of the adiabatic theorem and the pre- 
cise conditions for its validity. 

Despite its practical importance, no consensual and 
rigorous necessary and sufficient conditions for the va- 
lidity of the adiabatic theorem have been given. Only 
recently a proof that the commonly used textbook con- 
dition [l| is necessary for non-degenerate Hamiltonians 



[B| but not sufficient was given. For degenerate sys- 
tems, even a clear presentation of the theorem is lacking, 
let alone necessary and sufficient conditions. It is this 
paper's intention to fill that gap. 

With that in mind, our goal is three fold. First, us- 
ing techniques developed in @, Q , we aim at providing 
a clear and rigorous version of the adiabatic theorem for 
Hamiltonians with non-degenerate and degenerate spec- 
tra using a single formalism. We want to be as precise as 
possible in stating the adiabatic theorem to avoid com- 
mon misunderstandings Q , mainly due to a lack of quan- 
titative rigor in the way the theorem is usually presented. 
Second, we prove necessary and sufficient conditions for 
the validity of the rigorous version of the adiabatic the- 
orem here presented. The necessary condition for de- 
generate spectra reduces to the one in Q when no de- 
generacy is present. To obtain a sufficient condition, we 
rely on the adiabatic perturbation theory developed in 
[H Q. Finally, we apply these ideas to an exactly solv- 
able time-dependent degenerate problem where we 
show that the necessary and sufficient conditions here de- 
veloped provide the correct conditions under which the 
adiabatic theorem holds. 

To properly formulate the degenerate adiabatic the- 
orem (DAT) we first need to introduce the degenerate 
adiabatic approximation (DAA). As we will see, DAT is 
essentially a statement about the mathematical condi- 
tions for the validity of DAA. This understanding of the 
essence of the adiabatic theorem is akin to the ones of 
Berry [To[ and Tong [H[ , for non-degenerate systems, and 



111 and Wilczek hj, for 



to the ones of Wilczek and Zee 
degenerate systems. 

Degenerate Adiabatic Approximation. Consider an ex- 
plicitly time-dependent Hamiltonian H(i) with orthonor- 
mal eigenvectors \n 9 " (t)), where g n = 0, 1, • • • , d n — 1 la- 
bels states of the degenerate eigenspace rl n of dimension 
d n and eigenenergy E n (t), H(i)|n»»(t)) = E n {t)\n^(t)); 
and assume that d n does not change during the total time 
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evolution, t G [0,T]. An arbitrary state at t = can be 
written as |*(°)(0)) = £„ E^o &n(0)£^(0)K»(0)), 
where |6„(0)| 2 gives the probability of the system be- 
ing in cigenspace "H rl and \b n (0)U% n (0)\ 2 the prob- 
ability of measuring a specific eigenstate. A given 
initial condition within an eigenspace is characterized 
by one value of h n = 0, 1, • • • ,d„ — 1. A compact 
way of representing all possible initial conditions span- 
ning the orthonormal eigenspace H n is |\I> ( ' ' ) (0)) = 
E„=o& n (0)U»(0)|n(0)), where \n(t)) = (\n°(t)), \n\t)), 
. . . , \n dn 1 (t))) is a column vector, and U n (0) 
unitary matrix, U"(0)(U n (0)) t = 1. A particular initial 
state corresponds to choosing the corresponding element 
of the column vector |\l>^(0)). 

Then, the most general way of writing DAA is 
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(t))=^e-^W6„(0)U»(i)|n(t)), (1) 



n=0 



where oj n (t) = J Q E„(t')dt' /% is the dynamical phase, and 
the unitary matrix U n (t) = U n (0)T exp(/ * A nn (t')dt') 
the non-Abelian Wilczek-Zee phase (WZ phase). Here 
T denotes a time-ordered operator, and A™ n g (t) = 
(M£™ m (£))* a d n x d n matrix defined as 



[M mn (t)~ 



g m h„ 



(n h »(t)\m°™(t)), (2) 



with the dot meaning time derivative. For example, for 
a system starting at the ground eigenspace (b n (0) = 8 n o) 
|* (0) (i)> = e- iw °WU°(t)|0(t)). 

The time evolution of an informationally isolated quan- 
tum system is dictated by the Schrodinger equation (SE) 
ifi|*(t)) = H(i)|*(t)). What are the constraints on the 
rate of change of H(f) under which the system's evolved 
state \*&(t)) gets close to DAA? The adiabatic theorem 
we formulate next sets the conditions under which DAA 
holds. In other words, it precisely states when the sys- 
tem's dynamics can be approximated by DAA. 

Adiabatic Theorem: If a system's Hamiltonian H(t) 
changes slowly during the course of time, say from t = 
to t = T, and the system is prepared in an arbitrary su- 
perposition of eigenstates of H(t) at t = 0, say |\E^ )(0)), 
then the transitions between eigenspaccs T-L„ of H(<) dur- 
ing the interval t G [0,T] are negligible and the system 
evolves according to DAA. 

The three important concepts, slow, negligible, and 
evolved state, need further explanation. First, DAA is 
based on the assumption that the rate of change of H(i) 
is slow. A crucial matter is then to establish the mean- 
ing of slow precisely. Intuitively, the latter notion can 
be understood as a relation between a characteristic in- 
ternal time of the evolved system Ti, encoded in H(i), 
and the total evolution time T, such that Ti/T <C 1. For 
a fixed and finite Ti, one can always choose an evolu- 
tion time T that satisfies this condition. This state of 
affairs, however, is not satisfactory from a mathematical 



standpoint. Indeed, a main source of controversy in the 
literature arises from the lack of a precise quantification 
of the term slow. By using the degenerate adiabatic per- 
turbation theory (DAPT) Q, a generalization of APT 
Q , we can give a precise meaning to this notion of slow- 
ness, which is the key ingredient to the derivation of the 
sufficient condition of DAT. Second, to establish the nec- 
essary condition we follow Tong @ and others [IM3 
and assume that if the system's state is well described 
by DAA then all measurements performed on the system 
at any time must indeed be consistent with this assump- 
tion. This has a profound implication on the approximate 
dynamics the system obeys |5|. The following necessary 
and sufficient conditions provide the mathematical rigor 
required to make those concepts precise. 

The necessary condition. There is no unique way of 
establishing how close two quantum states are, imply- 
ing that there is no unique distance measure between 
states. A popular choice in the context of quantum in- 
formation is the fidelity measure. We stress though that 
DAT is not a statement about the fidelity between the 
true time-dependent state |*(<)> and DAA |* (0) (t)) be- 
ing close to one, i.e., \{^{t)\^^(t))\ ~ I. It is more than 
that, it is a statement about DAA expectation value of 
any observable being close to the exact ones. This notion 
is crucial to define geometric phases, thus for particle ex- 
change statistics, and is crucial for the philosophy behind 
DAPT and the proof of necessity that now follows. 

If DAA is an accurate description of the time evolution 
of a degenerate system starting, with no loss of general- 
ity, in its ground eigenspace (b n (0) = S n0 ) then = 
|* (0 V)) + 0(1/T) » |* (0) (t)), with \\0(1/T)\\ max « I, 
where 1 1 • 1 1 max is the max norm (the absolute value of the 
greatest element of a given vector/matrix). It immedi- 
ately follows that the system (a) approximately satisfies 
SE i7i|* 0)) w H(t)|*°(i)) which implies @ |*(f)) 
ISj^ ) (£)); and that (b) transitions to excited eigenspaces 
are negligible 0, || (n(t)| T |*(t)) T || max < l,n^ 0. 

Now, using (a), (b), and defining A nm (t) = E n (t) — 
E m {£) we notice that for n ^ [3] , 



<n(i)n*(t)> 5 



(n(t)\ T (H(t)~E (t))\*(t)) T 
A n0 (t) 

(n(t)\ T (ih\*(t)) T -E (t)\*(t)) T ) 



A n0 (t) 



ih(n(t)\ T \^°\t)) 



= ifar itJn(t) 



A n0 (t) 

n(i)| T [U°(t)|6(i))f 



A„ (t) 



where (n(i)| T |0(i)) T = 0. Taking the max norm on 
both sides and using (b) we get the necessary condition 
^||(n(t)r[U°(t)|6(t))] T /A„ (t)|| max « 1, n ± 0, t € 
[0,T\. Finally, using that ||U n (£)|| max < I leads to a 
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stronger WZ phase-free necessary condition, 



M n0 (t) 



Ano(i) 



<1, n^O, te[0,T], 



(3) 



where = max 5Zf =1 |fflij| for a, p x q dimensional 

i<j<p 

matrix A. When the spectrum is non-degenerate {d n = 
1), Eq. © reduces to the necessary condition of Ref. 

The sufficient condition. The first stcpt to stablish 
the sufficient condition is to prove the convergence of 
DAPT in its full generality. Intrinsic to the formulation 
of DAPT is a Taylor series expansion in terms of the 
parameter v = 1/T, and a necessary rescaling of time 
according to s = vt with s f= [0, 1] For small enough 
v one can always make DAPT converge (cf. Eq. ([6])). 

Inserting the ansatz 



(4) 



n=0p=0 



into SE with Cl p) (s) = e-i u "Wv p B% ) (s) and b£ p) (s) = 
Em=o e " U ° m ' S ^'»' l ( s )' DAPT gives recursive equations 
for Bto]j(s) in terms of lower order in p coefficients Q. 
The zeroth order is exactly DAA, with WZ phase nat- 
urally appearing as a requirement for the consistency of 
the series expansion. Note that for each n we have a se- 
ries involving the matrix Cn ( s )> P = 0, 1, . . . ,oo. The 
matrix element [C„ (s)]h n g n is the coefficient giving the 
contribution to order p of the state \n 9n (s)) to the solu- 
tion to SE. Here h n handles different initial conditions 
and for definitcness we pick the case h n = 0, Vn. Apply- 
ing the ratio test for series expansions, if the condition 



lim 

p— >oo 



[Cl P+1) (s)]o s „/[C^( S )] 0s „ <1, Vn,g ni (5) 



is satisfied for all coefficients then we guarantee conver- 
gence of DAPT. We can simplify further ([5]) by invoking 
the comparison test 



V Em=0 



lim 

p— >oc 



Em=0 



0g„ 



B (p) Cs") 



<>y„ 



Imposing that |[Cl P+1) ( S )] 0s J « \[CX' (s)] 0g J, 

\/n,g n , meaning that the zeroth order dominates, is 
equivalent to 



< 1, Vn,flVi 



i(0). 



(6) 
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p— m— 



p+l 



run V- 



«E [ B ^( 



0ff„ 



(7) 

which together with Eq. ((6|) are the rigorous sufficient 
conditions for the validity of DAA. In practice it is ex- 
tremely difficult to compute the previous limit when 
p — > oo and all orders p. We can come up, nevertheless, 



with some practical condition of convergence by looking 
at the ratio for a couple of finite orders p. Working with 
increasing p we get more and more conditions that, in 
the non -deg enerate case, can become stronger than the 
ones in |15j . In its simplest form, we may consider only 
p = 0. In this case both expressions merge into one and 
we demand it to be much smaller than the smallest non- 
null term appearing in the rhs of (|7jl . Thus, the practical 
sufficient test reads 



«E [ B ™«( s ) 



Off* 



<C mm 

n,g„ 



E 



0g„ 



(8) 



Using |<| B^l(s) = b n (0)U n (s)5 mn and the fact that 
at t = the initial state is |0°(0)) (6„(0) = <5„o) we get 



vj2 [b« oo 



m=Q 



Off, 



<C mm 

90 



[u°W] 



Offo 



,Vn,g n , (9) 



which is our intuitive and practical sufficient condition. 

, with 



3 -^m(s) 



B (1) <VI 



0g„ 



Indeed, noting that v E m =o 1 

ti ^ 0, gives the first order contribution of the excited 
state \n 9n (s)) to the wave equation, and that for n = 
it is related to the first order correction to the WZ phase 
0, it is clear that they must be much smaller than the 
smallest coefficient appearing in the zeroth order if we 
want DAA to hold. 

Equation ([9]) also depends on U' l (s) because B„i(s) 
depends on U n (s). However, a similar calculation to the 
one done for the necessary condition gets rid of these 
unitary matrices leading to [l4[ 



^(t)«min( [U°(t)] 



'0g 



t€[0,T], (10) 



where for n = and Vgo we have D° go (t) equals to 

Ei:~u\[M° n (t')(M^(t'))% >l0 \\ 




|A „(i')l 



and for n ^ and Vg n , D" (t) is given by 



(11) 



\do-l 



-l,d„-l 



|A n0 ( 



J2\ [M° n (i)] fcoff „| +d n J2\ [M 0n (0)] koln \\. (12) 



U'o=0 



fe ,;„=o 



Example. We now apply the previous ideas 
to a doubly degenerate four-level system subjected 
to a rotating magnetic field of constant magnitude 
B(f) = Br(t) and in spherical coordinates r(t) — 
(sin 9 cos wt, sin # sin if t, cos 6), with w > and < 
8 < ir being the polar angle. The Hamiltonian de- 
scribing this system is @, EH H(t) = Hbr(t) ■ T/2, 
where b > is proportional to the coupling between 
the field and the system and T = (T x ,T yi T z ) are the 
Dirac matrices Tj = a x ® o~j, j = x,y,z. Here Oj 
are the standard Pauli matrices implying the following 
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algebra for {Ti,Tj} = 2%J 4 , [T^Tj] = 2ie ijk U k , 
where J4 is the identity matrix of dimension four, Sij 
the Kronecker delta, the Levi-Civita symbol, and 
lift = I2 (8 Cfe. Starting at the ground state |0°(0)) the 
time-dependent solution in terms of the snapshot eigen- 
states is @ = e™*/ 2 [(1 + cos0) A_(t) + (1 - 

cos0)A+(i)]/2 |0°(t)) + e-™*/ 2 sin - A-(t))/2 

lO 1 ^)) + e hut/2 sin 2 (5+(t) + B_(t))/2 |l°(t)> + e^™*/ 2 
sin [(l+cos0)B-(t) - (l-cos0)5+(t)]/2 \l x {t)), wher e 
A±(i) = cos(ft±/j/2) + i (&±wcos0) sin(fi±t/2)/fi±, 
B±(i) = i w sm(fi ± f/2)/Q±, ft^ = w 2 + 6 2 ± 2wbcos6. 

Necessary condition. Since in this example Eq. 
© is [M 10 (i)] n = -[M 10 (i)] 00 = iwsin 2 (0)/2 and 
[M w {t)] w = -[M w {t)}^ = -iwsm(20)e iwt /4 the neces- 
sary condition ^ becomes w sin | sin0+cos0|/(26) <C 1. 
Our task now is to look at the exact solution, impose that 
DAA holds, and see if it implies the necessary condition 
above. If DAA holds then the absolute values of the co- 
efficients multiplying 1 1° (t ) ) and must be negligi- 
ble. This leads to Q wsm0f(6)/{2b) < 1, with f{9) = 
\b/VL+ + 6/0_ + cos0(6/O_ - b/Cl + )\. Noting t hat /(0) 
has a global minimum at 8 = ir/2 equal to 2b/\/b 2 + w 2 
it is not difficult to see that if w < b than /(0) > 
y/2. Hence, 1 > tusin0/(0)/(26) > w sin 0^2/ (2&) > 
wsin0| sin0 + cos0|/(26), which is exactly the necessary 
condition. When w > b we have /(0) > by/2/w which 
leads to wsin0/(0)/(26) > V2sin(0)/2 w sin0. Since 
sin0 « 1 DAA is not a faithful approximation to the 
exact state for general when w > b. This is expected 
since the rotating frequency w of the magnetic field must 
be much smaller than the coupling constant b (natural 
frequency of the system) for DAA to hold. The patho- 
logical situation where sin — > and the fidelity of the 
state approaches unity even though w > b does not lead 
to a state evolving according to DAA Q ■ 

Sufficient condition. Equations (|10til2[) become 
for t e [0,T], u. 2 £sin 2 (0)/6 <C min(|[U (i)] 09o |) 

90 

and 5wsin0(|cos0| + sin0)/(26) < min(|[U°(i)]o flo |), 

So 

where |[U°(i)]oo| = (1 - sin 2 0sin 2 (u;tcos(0)/2)) 1 / 2 and 
|[U°(i)] i| = sin0|sin(wtcos(0)/2)|, with g = 0,1. 
Note that the sufficient condition here is stronger than 
the necessary one because 5w sin 0(| cos 0| + sin 0)/ (2b) 
> u>sin0|cos0 + sin0|/(26). Moreover, looking at 
Eqs. §5§ and ©, and in particular (fT2"j) . we can 
show that in general the practical sufficient condi- 
tion implies the necessary one whenever the gap is 
constant. Since for this example the natural choice 
for the perturbative parameter v is the rotating fre- 
quency of the field (v = w) [IJ, we have wt < 1 for 
t g [0,T]. This implies that |[U°(t)]oo| > |[U°(*)]oi| 
and 5wsin0(|cos0| + sin0)/(26) > w 2 isin 2 (0)/6 dur- 
ing the whole evolution of the state. Hence, the 
sufficient condition boils down to only one equation, 
5wsin0(| cos0| + sin0)/(26) < sin0| sin(wt cos(0)/2)|, 
leading to 5w/(2b) < | sm(wt cos(0)/2)|/(| cos0| +sin0). 



Note that when t « Oand/or0 m tt/2, |[U°(t)] i| ~ and 
we must work with the non-null coefficient |[U (t)]oo|- In 
this case the sufficient condition is 5w/(2b) <C 1. 

It is important to remark now that if w > b we cannot 
satisfy the sufficient condition, no matter what the value 
of sin0 is. Indeed, since both terms appearing at the rhs 
of the sufficient conditions are smaller than one, assuming 
w > b leads to a lhs greater than one. The sufficient 
conditions are then consistent with the cases where the 
necessary condition fails. We cannot have sin0 w and 
w > b as an instance in which DAA holds. 

Our last task is to show that for w < b these condi- 
tions imply DAA. In other words, we must use them to 
show that the absolute values of the coefficients multi- 
plying |l°(t)) and of the exact solution are neg- 
ligible. Working with the largest of those this is equiv- 
alent to showing that [H| wg(9)/b <C 1, with g(9) = 
sin (6/f2+ + 6/0_) . Using that g(9) has a maximum, 
for e [0,7r], at = tt/2 given by 2b/(b 2 + w 2 ) 1/2 we get 
wg(9)/b < 2w/(b 2 + w 2 ) 1 / 2 < 2w/b. Hence, if the suffi- 
cient conditions imply that 2w/b <C 1 we are done. But 
noting that | sin(u)/jcos(0)/2)|/(| cos0| +sin0) < 1/2 the 
sufficient conditions reduce to bw/b <C 1 which obviously 
implies 2w/b <C 1. 

Summary. We established one rigorous necessary con- 
dition and two sufficient conditions, one rigorous and one 
practical, for the validity of the quantum adiabatic theo- 
rem for systems with degenerate spectra. Concepts such 
as "slowly/adiabatically changing Hamiltonians" and the 
"adiabatic approximation" for degenerate systems, of 
greatest importance for the implementation of adiabatic 
and topological quantum computation as well as non- 
Abelian fractional statistics, were quantitatively stated. 
It is this quantitative specification that allows for a pre- 
cise and rigorous formulation of the adiabatic theorem. 
Finally, we applied the adiabatic theorem to an exactly 
solvable degenerate problem, and provided a complete 
characterization of the mathematical conditions under 
which the degenerate adiabatic approximation holds. 
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